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ABSTRACT 

> 

r^ ' Context. LkHa 233 is a Herbig Ae/Be star with a collimated bipolar jet. As such, it 

^*i ' rnay be a high-mass analogue to the classical T Tauri stars and their outflows. 

t"*^ ■ Aims. We investigate optical forbidden lines along the LkHa 233 jet to determine physi- 

^h . cal parameters of this jet (electron density ric, hydrogen ionisation fraction x^, electron 

^^— ^ , temperature Tc). The knowledge of these parameters allows us a direct comparison of 

(— ^ , a jet from a Herbig star with those from T Tauri stars. 

Methods. We present the resuhs of HST/STIS and WFPC2 observations of LkHa 233 
(s^ " and its environment. These are the first observations of this object with a spatial res- 

$_( ' olution of < O'.'l at optical wavelengths. Our STIS data provide spectroscopic maps 

that allow us to reconstruct high angular resolution images of the bipolar jet from 
LkHa 233 covering the first ~ 2000 AU from the star in the blueshifted outflow lobe 
and « 4000 AU in the redshifted lobe. These maps are analysed with a diagnostic code 
that yields ric, Xc, Tc, and mass density nn within the jet. 

Results. The WFPC2 images in broad-band filters clearly show a dark lane caused either 
by a circumstellar disk or a dust torus. The circumstellar environment of LkHa 233 
can be interpreted as a conical cavity that was cleared by a bipolar jet. In this interpre- 
tation, the maximum of the optical and near-infrared brightness distribution does not 
coincide with the star itself which is, in fact, deeply extincted. In the blueshifted lobe, 
rie is close to or above the critical density for [SII] lines (2.5 x IC'cm"^) in the first 
arcsecond and decreases with distance from the source. The ionisation x^ ~ 0.2 — 0.6 
gently rises for the first 500 AU of the flow and shows two re-ionisation events fur- 
ther away from the origin. The electron temperature To varies along the flow between 
lO^K and 3 x 10* K. The nn is between 3 x 10^ and lO'^cm"^, and the mass flux 
M ^ 10"* - 10"^MQyr"\ The (radial) outflow velocities are si 80 - 160kms"\ and 
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they appear to increase with distance from the source. In the redshifted lobe, the ex- 
citation conditions are quite different: To, rie, x^, and nu are all lower than in the 
blueshifted lobe, but have the same order of magnitude. 

Conclusions. All these derived parameters are just beyond or at the upper limits of 
those observed for classical T Tauri star jets. This may indicate that the flows from 
the higher mass Herbig stars are indeed scaled-up examples of the same phenomenon 
as in T Tauri stars. 

Key words. ISM: jets and outflows - Stars: circumstellar matter - Stars: emission-line. 
Be - Stars: individual: LkHa 233 



1. Introduction 

Herbig Ae/Be stars (HAeBes) are intermediate-mass young stars, covering the mass range 
between 2 and 10 Mq. They are less well-studied than their lower mass counterparts, the 
T Tauri stars, for several reasons. First, they are rare because the initial mass function 
declines steeply with rising mass. Second, they spend much less time in their pre-main 
sequence phase than low-mass young stars. Third, they are embedded within envelopes, 
which makes observations of their immediate environment difficult. 

An interesting question is whether the environment of HAeBes is similar to that of 
lower mass young stars, for which the pre-main sequence phases are relatively well known. 
In particular one would like to know if ma ss accretion from circumstellar disks and mass 
ejection by me ans of jets and outflows (e. g 



Hartigan et al 



Ferreira 



1995; 



Konigl &: Pudrita . 



2000 



20021 ) is also at play in the formation of these higher mass stars. This would help 



in understanding if the accretion/ejection engine, a key element regulating the formation of 
low-mass stars, can be generalised to different stellar masses. Ver y high-mass young stars 



are known to possess disks and outflow s (see Zhang et al 



from brown dwarfs have been detected (jWhelan et al 



2005 



20051). and recently, also jets 



20071 ). To unveil the nature 



of the circumstellar surroundings of intermediate mass stars like HAeBes would thus allow 
us to add an important piece to the star formation puzzle. 



Indeed, some parsec-scale outflows from HAeBes are known ([McGroartv et al.l . 120041 ) 
similar to those from low-mass stars. To obtain information about the accretion/ejection 
engine, however, high angular resolution observations are necessary. The physical processes 
involved are confined to a region of at most 500 — 1000 AU from the central star. Therefore, 
even for the nearest star formation regions, sub-arcsecond resolution is necessary for ob- 
servations aimed at testing the accretion/ejection models. Following these ideas, we have 
observed the HAeBe star LkHa 233 and its environment with the Hubble Space Telescope 
(HST), which reaches an angular resolution of O'.'l in the optical. 

LkHa 233 has been classified as a HAeBe (spect ral type A7e) in the orig inal paper by 



Hernandez et al. 



(120041 ). the star has 



Herbig] (jl960f ) that defined the group. According to 

a mass of M^, ~ 2.1 Mq, a luminosity of L^, = 21.9Lq, and an age of about 7 Myr. These 



Send offprint requests to: Stanislav Melnikov, e-mail: melnikov@tls-tautenburg.de 

* Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the 
Space Telescope Science Institute, which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS5-26555. 



Stanislav Melnikov et al.: A HST study of the Herbig Ae/Be star LkHa 233 



parameters were calculated with a normal extincion law (i. e. R = 3.1 and A 



Y_ 



2^3). 

The extinction, however, is quite ambiguous in this region (jHernandez et al.l . 120041 ) and 
real parameters may be different from those. This star is located within a dark cloud that 
contains several other young stars. Throughout this paper we will assume the distance 
of LkHa 233 to be 880 pc and its heliocentric radial velocity to be —17.2 kms~^. These 
quantities were m easured for HP 213976 , which i s a bri ght B1.5V star associated with the 



same dark cloud (Calvet fc Cohen 



1978; 



Wilson 



19531 ). LkHa 233 is embedded within a 
reflection nebula that shows spikes at position angles PA = 50°/230° and 90°/270°. A 
first high ang ular resolution study of the central source and its close environment has been 

(jl993r ). Using speckle interferometry in the near- infrared, they 



Leinert et al. 



presented by 

obtained inform ation on a sub-arcsecon d scale and resolved a halo that has a radius of 

about 1000 AU. ICorcoran &: Ravi (|l997l ) noticed the presence of strong optical forbidden 



line emission from LkHa 233, indicative of a small-scale jet from this star. This jet was 
indee d discovered later by the same authors using longslit spectroscopy (jCorcoran fc Rav , 
19981 ). We observed the small-scale bipolar jet associated with this source with the Space 
Telescope Imaging Spectrograph (STIS) and have taken additional data obtained with the 
Wide Field Planetary Camera (WFPC2) from the HST Archive. 

The analysis of these data has led us to determine the structure of the jet at its base 
on small scales, allowing a comparison with analogous outflows from T Tauri stars in the 
same region of the plasma jet. In addition, the applicati on of a well-tested spectros copic 
technique for jet diagnostics to the observed line ratios (jBacciotti fc Eisloffell 119991 ) has 
provided values for key physical quantities along the jet at high spatial resolution, a step 
that is necessary for identifying the dynamics of the outflow. The obtained results point 
toward a strong similarity between the accretion/ejection engine in T Tauri stars and 
HAeBes, opening interesting prospects for a generalisation of the mechanism to all YSOs 
(Young Stellar Objects). 

In our paper, the methods of observations and data analysis are presented in Sect. [51 
and the results are illustrated in Sect. [31 In Sect. [H we discuss the physical implications of 
our findings and compare our results with previous studies of jets powered by lower mass 
YSOs. Finally, in Sect. [5l we summarise our conclusions. 



2. Observations and Data Analysis 

The jet of LkHa 233 was observed with the STIS spectrograph on board the HST on 
4 Oct 1999. The method of data acquisition was virtually identi cal to HST/STIS ob - 
servations of t he small-scale jets fr om the T Tauri stars DC Tau (JBacciotti et al.l . 120001 ) 



and RW Aur (Woitas et al 



20021) . The spectrograph slit was moved across the jet in 



seven st eps (S1...S7) of 0"07 from NW to SE, keeping it parallel to the jet axis (PA 



248° 



McGroartv et al.ll2004l ). The observations made use of the G750M grating that. 



in addition to Ha, covers the most prominent forbidden emission lines [01] AA6300, 6363, 
[Nil] AA6548, 6583, and [SII]AA6716, 6731. The pixel scales are 0.554 A/pixel and O'.'05/pixel 
for the dispersion and spatial direction, respectively. To ensure optimum positioning, the 
slit was peaked on the source before offsetting for each observation. As we show in the 
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following, however, the maximum of the optical brightness distribution does not coincide 
with the star. As a consequence, the central slit position, 5*4, does not coincide with the jet 
axis. The jet axis is instead closer to 53, and 5*6 and 5*7 unfortunately contain no signal 
from the jet. 

The basic data reduction steps were carried out using the HST/STIS pipeline. Cosmics 
and bad pixels were removed from the frames. To see the jet in the forbidden emission lines, 
the continuum contribution from the surrounding reflection nebula had to be subtracted 
carefully. This was done by fitting the continuum around the lines in each pixel row with 
a high-order polynomial and then subtracting it from the frames. The spectra were then 
used to reconstruct 2D images of the jet in the various lines and to diagnose the thermal 
conditions of the gas through our spectral diagnostics involving line ratios. 

The errors on the line fluxes were estimated considering the contribution of noise and 
uncertainties on the extinction. Ay. The noise was estimated from the subtracted back- 
ground as running averages of boxes of 7 x 8 pixels in size moved along both sides of the 
lines. The lines included in our STIS setting are too close in wavelength to allow a reliable 
measure ment of the extin c tion. An estimate in the region around this star is given, how- 

(|2004l ) who report a range between 2™3 and 3™ 7. To evaluate 



Hernandez et al. 



ever, by 

the possible effect of reddening, the lines were dereddened using the median value, and the 
uncertainty of Ay was estimated from the limiting values. The reddening uncertainty was 
added to the noise, but turned out to be only a few percent of the latter. Despite the fact 
that we could not measure the variation in extinction along the jet beam, one can therefore 
expect the resulting flux uncertainty to be negligible. Thus our further analysis was done 
using the measured flux values. 

HST/WFPC2 images of LkHa 233 were taken on 23 Jan 2001 in the broad-band filters 
F606W and F814W (Proposal ID 8216, P. I. Karl Stapclfeldt). After receiving the pipehne- 
reduced files from the HST Archive, cosmics and bad pixels were eliminated. Figure[l] shows 
images, with the target on the PC chip of the camera, with a scale of O'.'046/pixel. The 
achieved angular resolut ion is given b y the FWHM of the PSF that can be calculated using 
the Tiny Tim software (JKristl . Il995l ) and is ^ 0"! for the filters used. Previous adaptive 



optics observations of LkHa 233 presented by 



Perrin et al 



( 20041) have an FWHM of 0'.'27. 



3. Results 



3.1. The reflection nebula 

The WFPC2 images in both broad-band filters clearly show a bipolar ne bula with a dark 
lane b etween the two lobes (Fig.[T]). This lane is already indicated by the ICorcoran fc Rav 
(|l998f ) finding that the redshifted part of the bipolar jet is not visible for the first 0'.'7 
from its origin, wh ich is probably due to circumstellar dust occulting the receding flow. 



Perrin et al 



( 20041 ) have observed the dark lane on a larger spatial scale (FWHM = 0'.'27) 
using polarimetry in the near- infrared. The bipolar structure seen in Fig.[T] (left panel), 
together with the 'spikes' mentioned above, ca n thus be interpret ed as a conical cavity 
cleared by the jet, as has also been suggested by 



Perrin et al. 



(|2004l ). The optical emission 



comes from the inner edge of this cone that is illuminated by the star. An important conse- 
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Fig. 1. HST/WFPC2 image of LkHa 233 in the broad-band filter W814 (left panel). The 
axis of the bipolar jet (PA = 248° for the blueshifted lobe) is indicated. The jet is not visible 
in this image because the signal is dominated by continuum emission from the reflection 
nebula surrounding the star. The maximum of the brightness distribution whose size scale 
is 0'.'5 X O'.'S (marked with the white arrow in the left panel) is magnified in the right panel. 



quence of this is that the optical and near-infrared brightness maximum does probably not 
coincide with the star itself, but with a region on the surface of the cone pointing towards 
us. This is supported by the finding that this maximum appears elongated and does not re- 
semble a stellar PSF (Fig.[Tl right panel). The interpretation of the bri ghtness maximum as 
a scattering surface can also explain the extended structure found bv lLeinert et al.l (jl993l ) 
that was interpreted as a halo around the star in that work. From the present data, it is 
not clear if the 'dark lane' is a circumstellar disk seen close to edge-on or a dust torus 
around the star. If the former interpretation is true, LkHa 233 is one of the few young 
stellar objects where a disk has been directly observed. This disk, however, would have an 
apparent radius of about 1000 AU (adopting a size of the lane of approximately I'.'l at a 
distance to LkHa 233 of 880 pc, see Fig.[Tl left panel) and so it would be much larger than 
disks around classical T Tauri stars. The location of the bipolar jet is indicated in Fig.[TJ 
It is, however, not visible there because the continuum emission from the reflection nebula 
overwhelms the line emission from the jet. 



3.2. Jet structure 

Spectroscopic data obtained with HST/STIS in the way described in Sect. [H allow in prin- 
ciple to construct high angular resolution images of the flow in distinct velocity intervals. 
This is done by selecting the same velocity bin for each line in the S1...S7 spectra and ag- 
gregating the selected seven pixel columns in a 2D image. The result is analogous to what 
would be obtained with a high angular resolution integral field unit. This has been done 
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Fig. 2. Images of the LkHa 233 jet in different emission hnes, reconstructed from 
HST/STIS spectra, integrated over the full radial velocity range of the lines. The thick 
solid line marks the peak of the optical brightness distribution in the reflection nebula 
surrounding the star. The redshifted outflow lobe (PA = 68°) is on the right side of the 
image. Each panel has a spatial width of « O'.'S. The intensity is displayed logarithmically 
from 5 X 10~^^ to 3 x lO^^"* ergs^^ cm^^ arcsec^^. For Ha in the redshifted lobe, the signal 
is heavily affected by contributions from the disk/envelope in the first O'.'S from the central 
position. The bottom panel shows the total flux distribution coadding all forbidden lines 
along the jet. 



successfully in t he pa st for the jets fro m the classica l T Ta uri stars (CTTSs) DG Tau by 



Bacciotti et al 



pood ) and RW Aur by 



Woitas et al 



( 20021 ). In the present case, however, 



the signal-to-noise ratio of the individual spectra of LkHa 233 is rather low. Therefore we 
could not reach a similar level of detail for this jet, but we aimed at constructing high 
angular resolution maps integrated spectrally, by coadding (after subtraction of the back- 
ground) the column pixel values in each spectral channel. The resulting fluxes were strong 
enough to produce reconstructed images in each of the emission lines integrated over the 
whole radial velocity range covered by the lines with a sampling of O'.'05/pixel. Note that, 
since the STIS slit was moved from NW to SE, the numbering of the slit positions goes 
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from SI at the bottom pixel row of each panel to S7 at the top pixel row. The images 
in the emission lines included in the spectra are shown in Fig. [21 The lines are integrated 
over a velocity range of « [0, 225] kms^^ for the redshifted lobe and w [0, —225] kms^^ for 
the blueshifted lobe, and the intensity varies logarithmically from 5 x 10^^^ to 3 x lO^^"* 
ergs""'^ cm~^ arcsec"^. To make positions of the jet knots more evident, we show in the 
bottom panel the total flux distrubution coadding all forbidden lines along the jet. 

The redshifted lobe of the jet (PA = 68°), displayed in the right part of the image, 
is longer, but considerably fainter than the blueshifted outflow that occupies the lower 
part. The thick solid line marks the peak of the optical brightness distribution, which, 
as mentioned above, represents a bright point of the reflection nebula, and it is slightly 
displaced from the position of the star. This also caused the jet not to be centred on each 
individual panel, as explained above. 

The blueshifted jet lobe is clearly visible in the STIS data within its flrst two arcseconds. 
Five jet knots (labelled A...E in Figs.mandO can be resolved in this region close to the star, 
with possibly an additional faint knot (that we labelled Da) between D and E. The knots 
are seen more clearly in [01] A6000, Ha, and [SII] A6731, but they can also be distinguished 
in [Nil] A6583, where they appear to be more concentrated toward the jet axis. This is 
expected, since the [Nil] lines are a tracer of higher excitation which, in turn, is expected 
to be found in the axial region, where the jet shocks are more powerful. If an arbitrary jet 
inclination of 45° is assumed, the knot separations from the presumed position of the jet 
source are 200, 600, 1100, 1300, and 2100 AU, respectively. There is general consensus on 
the fact that bright knots trace the radiative cooling regions behind shock fronts, which 
form in the beam as a consequence of pulsation of the ejection mechanism at the jet 
base and/or internal instabilities during the propagation. The distance of the jet from 
Earth, however, does not allow the STIS resolution to ascertain whether the knots are 
actually bow-shaped, as seen in some other YSO jets. For the redshifted counter-jet, four 
jet knots (labelled A'...D' in Fig.[l]) are visible in the first 4" of the flow in [01] A6300 and 
[SII] AA6716,6731. The S/N ratio is, however, quite low here. We note that the knots in the 
two lobes are not placed symmetrically with respect to the source position, as one would 
expect if, e.g. the knots were due to instabilities in the accretion/ejection events. Similar 
asymmetries are fou nd in many other bipolar jets, such as e.g., the flow from RW Aur 



(Woitas et al 



20021 ) ■ but an explanation for these asymmetries is still lacking. Regarding 



the map in the Ha line, we note that the emission of the jet was mixed with scattered light 
contributions from the disk and the envelope. We thus used the signal from 5*6 and 5*7, 
which contains no jet emission, to attempt a subtraction of these contributions. However, 
the result should be taken with care, especially within the first 0'.'5 of the redshifted lobe. 
Here the signal appears to be heavily affected by the contributions from the disk/envelope 
and may not represent structures within the jet. 



Th e same jet has been observed at much coarser spatial resolution by 



Corcoran fc Rav 



( 19981 ). who spectrally resolved several velocity components and determined corresponding 
electron densities, as shown in their Fig. 2. Taking the difference in spatial resolution into 
account, their results are in good agreement with our derivations. In Sect.|3]we will compare 
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Fig. 3. Preliminary analysis of the ratios of selected forbidden lines along the blueshifted 
lobe of LkHa 233, integrated across the full velocity range and across the jet transverse 
section. Distances are from the brightness maximum of the circumstellar reflection nebula, 
here referred to as 'source'. Top panel: [01] A6300/A6363 with its error, compared with the 
theoretical ratio of Einstein's A coefficients (dotted line). Bottom panel: same as above, for 
the [Nil] A6583/A6548 ratio. In both panels discrepancies begin to be more evident from 
1'.'2 from the source. 

the images shown in Fig. [2] with analogous data resulting from high angular resolution 
analyses of the base of CTTS jets from other sources. 



3.3. Physical quantities along the flow 

Thanks to YSO jets possessing a rich spectrum of emission lines, one can determine many 
interesting physical quantities associated with the emitting gas. Besides the accurate mea- 
surements of the radial velocity through Doppler shifts, one can rely upon recently devel- 
oped spectral diagnostic techn iques to find, for example, the ionisation fraction , the total 



density, and the jet mass flux (JBacciotti fc Eisldffell . 



1999 



Hartigan et al. 



1994. 



To find values of these q uantities in the LkHa 233 jet, we follow a pr ocedure desc r ibed i n 



Podio et al. 



( 200611 



Bacciotti fc Eisloffell (|1999I ). and subsequently refined and extended in 
usually referred to as the 'BE technique ', which uses a combination of the forbidden lines 
from oxygen, sulphur, and nitrogen. These lines, excited coUisionally, can be considered 
optically thin in the density and temperature regimes found in stellar jets and are easily 
modeled. The Ha emission, albeit very strong, is not used in the diagnostics, because 
it is often contaminated by the reflection nebulosity around the star, it can suffer from 
absorption in high density zones, and it is produced by both coUisional excitation and 
radiative recombin ation in regions not spatially coincident with the emission region of the 
forbidden lines (see iBacciotti fc Eisloffell Il999f) . 



3.3.1. Preliminary analysis of the line ratios 

Prior to the diagnostic analysis, the spectra were carefully examined, and it was found that 
the redshifted lobe of the jet is too faint in our images for the diagnostics to be applied 
at high angular resolution with a reasonable level of confldence. Thus for this lobe we 
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integrated the signal spectrally and spatially in each of the four most prominent knots, 
and we applied the diagnostics to each knot as a whole (see Table 1). 

For the brighter blueshifted lobe, the S/N was also not sufficient to perform a 2D 
diagnostic analysis. Therefore, we integrated spectrally and across the jet, and the high 
angular resolution information was only retained for the direction along the flow. For 
the blue-shifted lobe, we then formed various ratios of the forbidden lines to check for 
the presence of evident problems of S/N, that may render these lines unusable in certain 
locations. The ratios [01] A6300/A6363 and [Nil] A6583/A6548 should equal the ratio of 
Einstein's A coefficients for spontaneous emission, which is about 3 in both cases. Departure 
from this value can give an idea of the regions of the jet where the diagnostics will be less 
reliable due to poor S/N. The comparison is shown in Fig. [3] with the theoretical value, and 
the line fluxes were integrated across the jet and in radial velocity, with errors propagated 
quadratically. 

The comparison shows that the observed points are fairly close to the theoretical line 
until about 1'.'2 from the source, corresponding to the tip of knot D in Fig[21 Beyond 1'.'2 
and up to 1'.'7, corresponding to a faint region of the jet, the observed ratio deviates from 
the expected value. The ratio of oxygen lines is acceptable again at the end of the jet at 
the location of knot E, while for nitrogen lines the ratio is still far from the expected one at 
this position. Our test implies that the diagnostic results for the regions of the flow farther 
from the source than knot D should be taken with caution, especially for the ionisation 
fraction, which is tightly correlated with the [Nil] measurements. Since it is likely that 
these problems are caused by the very low S/N of the faintest lines in each of these pairs, 
we decided to use the flux of the brightest one and the theoretical ratio for our subsequent 
excitation analysis. 

The top panel of Fig. [H illustrates the values of the ratio [SIIJA6716/A6731. From this 
ratio one can measure the electron density along the flow. The theoretical value of the ratio 
calculated for To = 10* K in both the low and high density limits are superimposed on the 
datapoints. The critical density for this ratio (corresponding to the high density limit) is 
2.5 X 10* cm~'^, while ric = 50 cm~"^ is the minimum value for the subsequent diagnostics 
to be applicable. From the plot, one can argue that the electron density along this jet is 
very high along the first 1" from the source, while it clearly drops at greater distances, 
despite the larger measurement errors. Where the measured [SII] A6716/A6731 values are 
outside the indicated range (because of poor S/N), the data can only provide upper/lower 
limits to the electron density, indicated in the following plots with arrows. 

In Fig.m we also report the observed value of the ratios 
[NII]A(6548+6583)/[OI]A(6300+6363) (hereafter referred to as [NII]/[OI], middle 
panel) and [01] A(6300+6363)/[SII] A(6716+6731) (hereafter [OI]/[SII], bottom panel), 
integrated in velocity and across th e jet cr oss section, with errors propagated quadratically. 
As shown in iBacciotti fc Eisloffell (jl999l ). the [Nil]/ [01] ratio only depends weakly on 



temperature, but it is a good indicator of the level of global ionisation in the flow, since 
the ionisation fractions of N and O are tightly correlated to that of hydrogen through 
charge exchange reactions. The plot shows a gradual increase in this ratio, and hence of 
ionisation, going from the source to the position of knot C. Then the values decline to a 
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Fig. 4. Preliminary analysis of the line ratios used for the diagnostics along the 
blueshifted lobe of LkHa 233. Fluxes are integrated across velocity and spatially across 
the jet transverse section. Distances are from the brightness maximum of the circum- 
stellar reflection nebula. Top panel: [SII] A6716/A6731 ratio with its measurement er- 
ror. The theoretical value in the low density limit (dotted line) and high density limit 
(dot-dashed line), calculated for To = 10'' K are superimposed. Middle panel: ra- 
tio [NII]A(6548+6583)/[OI]A(6300+6363), with its error. This ratio is very sensitive 
to the hydrogen ionisation fraction (see text). Bottom panel: variation of the ratio 
[01] A(6300-f 6363)/[SII] A(6716-f 6731), which is sensitive to both the electron temperature 
and the ionisation fraction of the emitting gas. 

minimum in the region of faint flux around 1'.'4 from the source. The ratio then appears 
to increase substantially at the end of the flow; however, this may not be real due to 
the aforementioned measurement problems for the [Nil] line fluxes. The [01] /[SII] ratio 
depends on both the hydrogen ionisation fraction and the electron temperature in the gas, 
but is correlated more tightly with the latter. The values of the ratio flrst decrease from 
the source to a minimum at 0'.'2, followed by a gradual increase along the jet, with the 
highest values reached at the position of knot C, and a minimum immediately downstream 
of knot D. The values of this ratio, however, are defined better toward the end of the flow 
with respect to the [Nil]/ [01] ratio. 



3.3.2. Results of spectral diagnostics 

In this paragraph we pres ent the results of the BE diagnostic procedure described in 
Bacciotti fc Eisloffell (|l999l ) applied to our dataset. Basically, the various quantities are 
obtained by comparing the observed line ratios with those calculated by a code designed 
to predict the fluxes in forbidden lines starting from a grid of values of electron density 
ric, hydrogen ionisation fraction x^, and electron temperature T^. The electron density is 
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obtained by inverting the ratio of sulphur hnes, which is independent of a;o and only weakly 
dependent on Tc. Then, Xc and To are found by combining and inverting the [01] /[Nil] and 
[OI]/[SII] ratios, also using the retrieved ric. The calculation was performed with a numer- 
ical code developed by one of us (U. Locatelli) to apply the technique to large datasets. 
For more details on the diagnostic techni que and on the u nderl ying physica l assumptions. 



together with its limits of application, see iBacciottil (|2002l ) and iPodio et al.l (|2006l ) 



Due to the aforementioned problems of poor S/N, we performed the diagnostic analysis 
on the ratios obtained integrating the spectral datacube both in velocity and in the spatial 
direction across the jet. Thus high angular resolution is retained only for the direction 
along the flow. We stress that, in the present case, due to the problems of low S/N, we 
only use the brightest lines of the [01] and [Nil] doublets to form the observed ratios that 
will be compared with the theoretical ones in the diagnostics (see Sect. I3.3.ip . 

The results of the diagnostics are collected in Fig. [5] The figure, in the first two panels, 
also illustrates the intensity profile of the emission lines and the radial velocity of the 
material at each position along the jet, to give a positional and kinematical reference for 
the values of the thermal properties shown in the other panels. 

The line fiux distribution along the flow (top panel) is derived from the sum of the 
spectroscopic images in [01] AA6000,6363, [Nil] AA6548,6583, and [SII] AA6716,6731, and 
subsequent integration of the combined image across the jet axis. As mentioned above, this 
step was necessary to reach a sufficient S/N ratio for application of the diagnostics. 

The information about the jet structure is complemented by determining of the radial 
velocity of each feature in the [SII] A6731 and the [01] A6300 lines (second panel from top). 
This was obtained by summing the spectra in the direction transverse to the jet, which is 
coadding the individual 5*1. ..57 and then performing Gaussian fits to the integrated (spa- 
tially) line profile in each pixel row. This yielded peak radial velocities versus distance from 
the source, which were then transformed to systemic velocities assuming the heliocentric 
radial velocity of LkHa 233 to be — 17.2kms^^ (see Sect.[T]). 

The velocities along the line of sight vary between —60 km s~^ close to the star and 
— 160 km s~^ at the end of the observed portion of the jet. Errors are not reported on each 
point to keep the plot readable, but these were accurately determined from the Gaussian 
fit and vary between 10 and 40 km s~^ for both lines, being greater toward the end of the 
flow, which is fainter. For the redshifted lobe, no precise Gaussian fit was possible due to 
the low S/N. A rough estimate of the line shifts, however, indicated velocities of the plasma 
of about the same amplitude as in the blue lobe. 

The velocity determinations indicate a gradual increase in velocity with distance in the 
examined jet portion. This increase is, however, not continuous, as the velocity is fairly 
constant along the first 0'.'9 of the jet (the brightest jet portion), and then presents sudden 
jumps at knots D and E. This probably indicates that the last two knots move into a 
more rarefied external medium. Internally to almost all of the knots, one can marginally 
distinguish a local decrease in the radial velocity across each knot, revealed in both lines. 
For an interpretation of this effect, see Sect.m 

The lower panels of Fig. [5] illustrate the thermal properties in the blucshifted lobe of 
the jet, resulting from the application of the BE technique. 
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Fig. 5. Physical properties of the emitting gas in the blueshifted lobe of the bipolar jet 
from LkHa 233 as a function of distance from the brightness maximum of the circumstellar 
reflection nebula. All quantities (except the radial velocity) are derived after integration 
across wavelength and across the jet width. From top to bottom; intensity profile resulting 
from the sum of forbidden line fluxes, in units of 10^^^ erg s^^ cm~^ arcsec^^ (letters refer 
to the jet knots displayed in FiglJ); radial velocity resulting from the analysis of the lines 
[01] A6300 and [SII] A6731, integrated across the jet (errors range between 10 and 40 km 
s~^, see text); electron density nc, hydrogen ionisation fraction Xc', electron temperature 
Te; total density tih = Uc/xc', mass loss rate Mj. 



The third panel from top describes the variation in the electron density, ric, along the 
jet, as derived from the inversion of the ratio of the two sulphur lines for an assumed 
temperature of 10^ K. The reported errors on tIc come mainly from the propagation of the 
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S/N uncertainty of the [SII] hnes, with a smaller contribution from the a priori uncertainty 
of Tc, although the line ratio is only weakly sensitive to this parameter. 'Upward' arrows 
are lower limits, equal to the critical density, in the positions where the ratio reached its 
saturation value. Similarly, 'downward' arrows indicate positions where the electron density 
is lower than 50 cm^^ (the limit of applicability of the BE technique). As already argued 
from the inspection of Fig.[?l the electron density is quite high in the first portion of the 
flow. Its value scatters around 10^ cm~^ in the first O'.'S and then gently decreases down to 
10^ cm"'^ at 0'.'75, corresponding to the base of knot C. It is then apparently extremely high 
inside this knot, even though the corresponding line brightness is not exceedingly large. 
The region downstream of knot C presents a different character; since it is much less dense 
in free electrons. Here n^ scatters between 500 to 1000 cm^'^, apart from those positions 
where the lower limit is reached, and an isolated point of higher density in knot E. Overall, 
the situation appears to be consistent with the jet experiencing less confinement by the 
surrounding medium downstream of knot C. 

In the fourth and fifth panels from top, we illustrate the results of the BE diagnostics for 
the hydrogen ionisation fraction Xc and the electron temperature Tc along the blueshifted 
jet. In determining these quantities, great attention has been devoted to a realistic estimate 
of the errors, as there are a number of different sources of uncertainty in this case. First of 
all, the diagnostic uses ratios of lines produced by different elements, therefore an assump- 
tion of the s e t of e lemental abundances is necessary. This aspect is extensively analysed in 
Podio et al.l (|2006r ). where the authors conclude that solar abundances are not adequate 



in all cases, but instead the abundance set appropriate for the region under study should 
be chosen. Unfortunately, no determination of abundances exists i n the literature for th e 



Lacerta Cloud, so we have adopted the abundance set of Orion by lEsteban et al.l (|2004l) 



assuming that the t wo regions have been exposed to the same enrichment processes (see 



Neufeld et al. 



20061 ). Another source of error comes from the uncertainty in rio, which is 



an input of the diagnostic calculation. We then have to consider the measurement error 
of the ratios [Nil]/ [01] and [SII]/ [01], as displayed in Fig.H which includes both the S/N 
accuracy and the effect of variations in extinction along the jet beam that are unaccounted 
for. 

All these sources of error were considered and their effect on the final determination 
of Xe and T^ was found a posteriori by running the code at the extremes of the range of 
variations in the input quantities. The resulting uncertainties are displayed as error bars 
in the panels. These are limited to 15-20% in the bright part of the jet, but they rise up 
to 50-60% on average in the faint final section. Missing points are due to the fact that in 
those positions the values of the ratios fall out of the range of validit y of the technique . 



This basically occurs in regions of poor S/N in the [Nil] lines, but see lPodio et al.l (J2006I ) 
for a more detailed discussion of this point. It has also to be kept in mind that the positions 
where the electron density was only determined as a lower limit are critical, and the results 
for Xo and Tc must be taken with caution in these regions. A better determination of the 
electron density in these positions would be desirable, but this would require the use of 
pairs of l ines of higher critical den sity not included in our STIS setting. To overcome this 
problem, iHartigan fc Morsel (|2007l ) have recently attempted a simultaneous inversion of all 



14 Stanislav Melnikov et al.; A HST study of the Herbig Ae/Be star LkHa 233 

the line ratios available from similar STIS spectra of HH 30 in the high density regions 
toward the star. The results for HH 30 are encouraging, but the mathematical complexity 
involved in implementing our code along the same lines does not seem justified by the very 
limited number of positions in this jet where the [SH] ratio gives poor information. 

Despite the many sources of error discussed above, our diagnostics allowed us to estab- 
lish with reasonable accuracy both the value and the variation in the ionisation fraction 
and the electron temperature in most of locations along the jet. In the first bright portion 
Xo appears on average to increase from about 0.2 in knot A up to 0.6 in knot C, where 
the jet shows the maximum ionisation. Note that the latter corresponds to a lower limit 
determination of the density, which may lead to a possible overestimate of the ionisation 
fraction. An inspection of Fig.[31 however, reveals that the [NII]/[OI] ratio itself presents a 
peak in that position, which lead us to believe that the enhancement of ionisation in this 
knot is real. Downstream of knot C Xe falls again to values around 0.3, but it clearly rises 
again in knot E up to 0.4 — 0.5. While the variations are moderate, they are nevertheless 
evident, and we note no clear association between peak intensity and peak ionisation inside 
each knot. 

The electron temperature To in the jet ranges between 1 and 3 x lO'' K, and its pattern 
of variation is complex, too. The value of To shows a marked decrease in the first 0'.'3 inside 
knot A, down to 10^ K at 0'.'25, and it then turns up again to 3 x 10** K at the end of 
knot B. The temperature shows a minimum in knot C, but here again the decrease to 
7000 — 8000 K may be artificially caused by the poor determination of the electron density 
in this knot. In knot D there are less positions in which the diagnostic is applicable. Here, 
Tc first appears to increase again to 3 — 4 x lO"' K at I'.'l from the source and then to 
decrease to 1.5 — 2 x lO'' K in knots Da and in the first part of knot E. There is a final 
local peak of the temperature to 2.2 x 10"* K toward the second half of knot E. 

Once Xc and ric are determined, a rough estimate of the average total hydrogen density 
riH can then be derived as nn = Uc/xc (errors are propagated quadratically from these 
quantities). Note, however, that this simple calculation likely overestimates the actual 
average densities in the jet by a factor 3 — 5, because the determination is biased toward 
the regions of higher line brightness, which in turn are those of higher compression after 
a shock front. Nevertheless, it is extremely useful to have at least an approximate value 
for this parameter, because the total density enters all the relationships governing the 
dynamics of the jet. In particular, it is extremely important to know this parameter when 
attempting a comparison between the analytical and numerical model predictions and the 
actual observations. The jet appears to be quite dense in the initial portion, as the hydrogen 
total density varies between 10'* and 10^ cni^-^ in knots A-D. Of course, the lower limits 
found for rto affect density estimations inside knot C. Downstream of knot D, jih presents 
lower values, scattering on average around 5 x 10"* cni^"^. Overall, the jet from this source 
appears to be one of the densest studied so far with the BE diagnostics. 

The last panel of Fig. [5] illustrates our estimate of the mass flux in the jet. This is 
discussed in detail in the next section. 

We conclude this part by commenting briefly on the diagnostic results obtained for the 
red lobe, from the data integrated spectrally and spatially in each knot (see Table [T|). The 
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Table 1. Physical quantities in the redshifted lobe of the jet from LkHa 233, averaged 
spatially and spectrally over each jet knots (see text). 



Knot Location log(nc) x^ log(nH) Mj 

(") (cni-3) (cm-3) (lO'^Mgyr-^) 



A' 


0.65 - LOO 


4.06 


0.12 


4.98 


3.7 


B' 


1.85-2.50 


3.85 


0.19 


4.58 


1.5 


C 


2.55-2.95 


3.94 


0.47 


4.26 


0.7 


D' 


3.00 - 3.45 


3.90 


0.20 


4.59 


1.5 



excitation conditions are quite different here from those found in the blue lobe, which is 
already indicated by the fact that this counter-jet is more prominent in [SII] AA6716,6731 
than in [01] A6300, in contrast to the blueshifted lobe. From the table one can see that ric, 
Xc, and njj arc all lower in the redshifted lobe than in the blueshifted one. In particular 
Tie < ncrit of the sulphur lines for all the knots investigated here. The maximum electron 
density is found in knot A', where it is about 10"* cm~^. It then decreases in the other 
knots, but it is never less than 8000 cm"'^. On the other hand, the ionisation fraction 
is generally not more than 10 — 20%, if exception is made for knot C, which seems quite 
excited, reaching Xc ^ 0.5. The combination of Xc and Uc provides quite a high total density 
that decreases from 10^ cm^'^ at the base of the jet down to ~ 1.8 x 10^ cni^"^ in knot 
C The electron temperature is around 1 — 2 x 10* K. Because of the low signal in [01] 
lines, however, the estimate is not as detailed as in the blue lobe, so values are not given 
in the table. We instead attempted an estimate of the mass flux, which is discussed in the 
following section. 

4. Discussion 

4.1. Mass flux in the jet 

One of the most useful quantities to know when comparing models with observations in 
the case of jets is the mass flux in the outflow, Mj, and its variation along the flow. For 
example, in the magnetohydrodynamic (MHD) models proposed to explain the formation of 
YSOs, the ratio between ejected and accreted matter around the star is fixed (Mj/Macc ~ 
0.01-^0.1). Moreover the knowledge of M,j allows us to estimate other important dynamical 
quantities like the linear and angular momentum carried by the jet (adding the information 
on poloidal and, where available, toroidal velocities). In turn, these estimates can clarify if 
the jet is capable of accelerati ng coaxial outflows , to clear the stellar environment and to 



inject turbulence in the cloud (jPodio et al.l . I200G ), as well as to ext ract the excess angular 



Woitas et al 



20051) 



momentum from the disk/star system (see e.g. 

In principle we can determine the mass flux from the emission lines in two ways, de 
scribed in detail in 



Nisini et al 



(|2005f ). In the first method, the mass fiux at each position 
is estimated by combining the total density inferred from diagnostic techniques with the 
values of the jet radius and velocities, rj and wj, provided by the FWHM of the luminosity 
profile across the jet and the kinematical analysis, respectively. This approach assumes 
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that each considered 'shce' of the jet is uniformly filled at the density derived from the 
diagnostics, which is biased toward more compressed regions. The second method com- 
pares the observed luminosities of the forbidden lines to the one calculated by an emission 
model, thus providing the mass of the emitting gas. This method implicitly takes the fill- 
ing factor into account, but is affected by uncertainties in absolute calibrations, extinction, 
abundance, and distance. In our case, the problems of an unknown filling factor is miti- 
gated by our observing the jet close to the source and at high angular resolution, which 
limits the possibility of including large regions of non-emitting material in the beam. On 
the other hand, extinction, elemental abundances, and distance are not well known in our 
case. Therefore, we prefer to use the first method, keeping in mind that the derived values 
of Mj have to be considered as upper limits. Previous studies performed on other jets by 
applying both methods show that the mass flux derived in this way is overest imated by 
not more than 20% when high angular resolution is used (JBacciotti et al.l . ll999f ). 

With our estimate of njj, the mass flux can be computed as Afj = /iTOpfiHTrrj wj, 
where rj is the jet radius, vj the full spatial velocity, rup the proton mass, and /i the mean 
molecular weight, equal to ~ 1.24 for interstellar atomic gas in standard composition. In 
our case, the jet radius was estimated taking half the FWHM obtained from Gaussian fits 
to the sum of all spectroscopic maps in forbidden emission lines, which is nearly constant 
at ~ 100 AU between d — O'.'l and 0'.'6. As the S/N is not high enough to derive reliable 
results beyond d — 0'.'6, we assume the jet radius to be the same there as closer to the origin, 
i.e. rj — 50 AU. The flow velocity vj is derived from fits to the coadded spectra S1...S7 
in [01] A6300 and [SII] A6731, as mentioned above, and de-projected with an arbitrary 
inclination angle of i = 45° . The resulting mass flux in each position of the blueshifted lobe 
is shown in the bottom panel of Fig. [5] (errors also include the uncertainties on vj and rj 
estimated from the Gaussian flt to the proflles of the emission lines) . The mass flux appears 
to be fairly constant at a value of about 10~^ Mq yr^^ between the source and knot G, 
indicating that the flow is well-coUimated in this portion and that turbulent interaction 
with the external medium, which may change the amount of transported mass, is negligible. 
The mass flux then decreases by 1 — 2 orders of magnitude in the last part of the jet. This 
could be because the real jet radius is actually larger than the one assumed. The fact 
that the jet undergoes lateral expansion downstream of knot C is also supported by the 
marked change of the other physical properties derived by the diagnostics. Alternatively, a 
significant portion of the mass could have been ejected sideways by the bow shocks in the 
working surfaces during the jet propagation, or Mj may be higher now than it was in the 
past. 

Regarding the redshifted lobe, for which our derivations come from data integrated 
spatially and spectrally on each con densation, w e can estimate the mass flux averaged 



over each knot as was done, e.g., in iPodio et al.l ()2006' ). To this purpose, we assume a 



de-projectcd flow velocity of 170 kms~^ and a jet diameter of 100 AU, similar to the 
blueshifted lobe. The obtained values are reported in Table 1. The resulting mass flux in 
the counter-jet seems to have the same order as in the blueshifted jet, with a maximum 
for knot A' , where we obtain Mj ^ 3.7 x 10~^Mq yr~^, and lower values in the following 
knots, but still the same order of magnitude. 
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Overall, the value of the observed mass flux in proximity to the star is slightly above 
the upper boundary of the values obtained for the jets from T Tauri stars. If the ratio 
of accreted-to-ejected mass flux is regulated as suggested by standard MHD jet launching 
models, the results would imply a mass accretion, Mace, onto the star of the order of several 



lO-^M, 



©yr" 



4.2. Nature of the knots 

In the second panel of Fig. [51 internal to almost all of the knots, one can marginally dis- 
tinguish a local decrease in the velocity, revealed in both the [01] and [SII] lines. This is 
consistent with the knots being mini- working surfaces, which should form inside the jet 
according to some models of jet propagation and are also directly imaged in resolution 
observations of T Tauri jets (JRav et al.l . 119961 ). Each working surface is contained within 
two shocks, one on the side of the jet source, the so-called Mach disk, in which the high 
velocity material of the jet beam coming from the source is decelerated upon entering the 
working surface. The other shock is on the downstream side, the bow-shock, in which the 
material ahead and moving more slowly than the shock front, is accelerated upon entering 
the working surface. The cooling region of the working surface is located between the two 
shocks and can be assumed to extend about 10^^ cm, which is the same scale as the knots 
observed in this flow. It is therefore reasonable to interpret the local variation in the veloc- 
ity observed within each knot as a consequence of the combined action of the two shocks 
constituting a working surface. In this framework, the slight differences in the position of 
the velocity variations in the two emission lines can be due to the peak emission region o f 
the two species behind a shock front being slightly displaced (JBacciotti fc Eisloffeill999f) . 
In conclusion, even if a bow-shaped morphology of the knots is not directly identified in 
the image of this flow, the velocity pattern appears to point toward a similarity of the knot 
properties with those of the TTS jets. 



4. 3. Comparison with flows from lower mass stars 

These findings can be compared with studies of other jets from T Tauri stars studied at 
high angular resolution and/or analysed with the same diagnostic method. The HH 30, 
DC Tau, and RW Aur jets were observed at high angular resolut i on with HST close to 



their sources (d < 1000 AUl (IRav et ah 



Woitas et al. 



2002 



adaptive optics techniques (JDougados et al 



1996; 



Bacciotti et al 



1998 



2000; 



Bacciotti 



2002; 



Hartigan fc Morsel. 12007) and from the ground with the CFHT using 



20001) . These jets were then analysed with the 



same diagnostic technique used here. Several other jets (HH 34, HH 1, HH 111, HH 46/47, 
HH 24, HH 83, HH 30, HL Tau, and Th 15-28) were observed from the gromid on much 
larger spatial scale s of 1000 to some 10000 AU from their parent s tars and ana l ysed with 



the BE technique (JBacciotti fc Eisloffell . 



1999 



Nisini et al. 



2005 



Podioetal 



2006h . In 



these cases the sources are either deeply embedded protostars or T Tauri stars. 

In the first w 2000 AU of the blueshifted LkHa 233 jet, the electron density Uc is 
remarkably higher than the one found in jets at large distance from the star, but similar to 
the electron density values found at the base of the jets observed at high angular resolution 
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( Bacciotti et al. 



2000 : 



Bacciotti . 



20021 ). On larger scales rio is expected to decrease due to 



recombination and the widening of the flow. 

In the HH 30 jet, analysed at high spatial resolution, the ionisation is lower on the whole 
than in LkHa 233. Likewise, it is rising in the first part of the flow and reaches its first 
maximum at d w 300 AU (in contrast to d w 500 AU for LkHa 233). In the same region the 
temperature decreases steeply from 2 x lO"* to less than lO"* K. Similar behaviour i s also seen 
in th e DG Tau and RW Aur jets observed with HST at high angular resolution (JBacciotti 



2OO2I Bacciotti et al. in preparation, Melnikov et al. in preparation), thus this appears 



to be a common phenomenon. Several studies have ind eed recently tried to justify these 
trend s in the framework of jet acceleration models (cf. 



Shang et al. 



2002 : 



Garcia et al. 



20011) . Alternatively, knot A might trace the presence of a strong shock, where To is high 



and Xc fi rst increases to reach a plat eau and then 'freezes' and slowly recombines along the 
flow (cf. iBacciotti fc Eisloffei Il999f ). For the other jets studied on larger angular scales, 
the ionisation is of the same order of magnitude as for LkHa 233. Re-ionisation events are 
found in a number of these outflows, but they occur on a much larger spatial scale than in 
the jet studied here. 

The total density in the LkHa 233 jet, derived by combining the ionisation fraction and 
electron density, is quite high, at least in its initial part, reaching 10^ cm~^. Overall, the 
jet from this source appears to be the densest studied so far with the BE diagnostics. 

The electron temperature To in the LkHa 233 jet is of the same order of magnitude as 
for the other jets. Peaks of To along the jet, as in the blueshifted lobe of the LkHa 233 jet, 
can also be seen in the outflows of HH 46/47, HH 24 C, and HL Tau, but again at much 
larger separations from their respective sources. Finally, the observed flow velocity is not 
very different from that of T Tauri jets. 

The mass loss rate is another fundamental parameter for testing the similarity of the 
physics in YSO outflows. In the case of LkHa 233 the mass loss rate is in general larger by 
up to one order of magnitude than for the jets observed on extended scales. It is, however, 
only slightly higher than the values found for other jets from lower mass stars studied at 
high angular resolution close to the source. In conclusion, the physical properties of this 
jet appear similar to those of classical T Tauri stars, apart from an enhancement of the 
tot al density. 



Garcia Lopez et al.l (|2006l) have determined accretion rates for a number of Herbig 



stars using a well-established correlation between T(Br^) and the accretion luminosity 
Lace- For a set of 36 Herbig stars, they found that 80% are accreting matter at rates of 
3 X 10"^ < Mace < IO^^Mq yr-i. If we adopt a mass flux for LkHa 233 of about 10% 
of its mass accretion rate (according to MHD models), then the Mace ^ IO^^Mq yr~^ 
(taking a smaller fraction of ejected to accreted matter would lead to an even higher value 
of Ma ce)- This would ma k e LkH a 233 one of the strongest accretors among the Herbig 
stars. iGarcia Lopez et al.l (J2006() also found that the Afacc distribution of Herbig stars is 
roughly consistent with the prediction of a Afacc oc Af^'* relation derived for subsolar mass 
objects and noticed a deficit of very strong accretors with rates higher than 10~^Mq yr~^. 
They suggested that this deficit could be due to a, i i "agi ng" -effect, since the accretion rate is 
expected to decrease with time (JHartmann et al.l . ll998l ). From this point of view LkHa 233 
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would be a very young Herbig star with a very high accretion rate in which the jet carries 
away an excess of angular momentum. 



5. Conclusions 

We have analysed images and spectra of the immediate environment of the Herbig Ae/Be 
star LkHa 233 taken in the visible range with the Hubble Space Telescope, which provides 
unprecedented high spatial resolution (O'.'l in the optical). The observed brightness distri- 
bution in archival HST images is probably from the surface of a cone; i. e. we do not see 
the source directly. We infer that a disk or dusty torus obscures the star. The brightness 
of the nebula does not allow one to see the jet in filtered narrowband (AA w 50 A) images, 
but it could be imaged in 2D maps derived from STIS data, reconstructing images of the 
flow in the various lines integrated in velocity from a set of seven spectra taken by stepping 
the slit in adjacent positions. 

In this way, we have resolved fine structure in the distant bipolar jet from LkHa 233 (880 
pc), in both the blueshifted and redshifted lobes. Within its first arcseconds (« 2000 AU), 
the jet knots appear circular or elliptical at HST resolution, although the kinematical 
analysis reveals that they probably trace mini-working surfaces, as is common in stellar 
jets. 

We have analysed t he spectroscopic maps with a diagnostic code designed to adapt the 



so-called BE technique (JBacciotti fc EislofFellll999r ) to large datasets. The low values of S/N 



of the spectra forced us to integrate the maps across the jet width in the blue-shifted lobe, 
and across the extension of each knot in the red lobe, in order to have acceptable inputs 
for the diagnostic technique. The results are displayed in Fig. [5] and Table [2 respectively. 
High angular resolution is retained along the jet axis in the blueshifted lobe. 

The electron density in the first arcsecond of the blueshifted lobe is about 10^ cm~^, 
close to the critical density for [SH] emission and, in a few positions, is above the critical 
value. The electron density decreases along the flow by slightly more than one order of 
magnitude, reaching values of about 800 — 1000 cm^'' at the position of the last visible 
knot (2" from the source). 

The hydrogen ionisation fraction in the same lobe varies between 0.2 and 0.6. It gently 
rises on the first 500 AU of the flow, like other jets from young stellar objects studied at 
high angular resolution; but instead of reaching a plateau followed by slow decay, it goes 
through apparent episodes of re-ionisation along the flow. 

The electron temperature T^, found with the diagnostics varies between 1 and 3 x 10** K, 
in agreement with the kind of lines observed. The jet knot A, located closest to the origin, 
is correlated with a maximum of the electron temperature and might represent a strong 
shock. In the region of jet knot C, there appears to be an anti-correlation between T^ and 
ionisation. The nature of this anti-correlation is unclear, but the result for the temperature 
here might have been affected by quenching of the [SH] lines. 

Combining the values of x^ and no we obtained estimates of the total density in this 
flow. The jet appears to be quite dense, especially close to the star; and with its tih ^ 10^ 
cm^"^, it turns out to be the densest outflow studied with the BE technique so far. 
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The mass flux Mj is « 10~^ Mq yr^^ in the first arcsecond of the jet, while it decays 
by an order of magnitude further out. These values are slightly higher than the mass-loss 
rate measured for the jets from lower mass YSOs. This is the first determination of a mass 
flux rate of a jet from a Herbig star obtained with the same method as for T Tauri jets. 

In the redshifted counter-jet, electron density, ionisation, and mass density are all lower 
than in the blueshifted jet and are consistent with a lower degree of excitation. 

We have compared the outcomes of our study with the properties of jets from T Tauri 
stars studied previously with the same diagnostics. The physical properties of the jet from 
this HAeBe star appear totally analogous to those in T Tauri jets, scaled toward higher 
densities and higher mass-loss rates. This appears to be consistent with a similar accre- 
tion/ejection engine in the two classes, with an enhanced rate of mass accretion onto 
HAeBes with respect to TTSs. 

Overall, these results suggest the same type of 'engine' drives the accretion/ejection in 
intermediate and low mass stars as well as sub-stellar objects. 
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